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We introduce and study the Wannier functions for an electron moving in a plane under the 
influence of a perpendicular uniform and constant magnetic field. The relevance for the Fractional 
Quantum Hall Effect is discussed; in particular it shown that an interesting Hartree-Fock state can 
be constructed in terms of Wannier functions. 
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I. INTRODUCTION 



It is widely accepted that Quantum Hall Effect can be described by a model of interacting electrons moving 
in a plane. For many features of the phenomenon the zero temperature limit seems to be a good approximation. 

In the attempt to describe the state in terms of single particle approximation and subsequent perturbations around 
this trial state, the author has introduced ||] a set of functions which are eigenfunctions of the single particle Hamil- 
tonian and moreover are invariant, up to a phase, under finite translations associated to a lattice. The phase depends 
on the coordinates. In subsequent papers iQ,^ an Ansatz was introduced in order to construct the trial state. In Ref. 
numerical evidence is given that this state yields a mean Coulomb energy close to the value provided by genuine 
' many-body states (e.g. Laughlin's 1^). From this property and from the fact that the trial state is given in analytic 
in . form in terms of a Slater determinant (or single particle approximation) the author hopes that a perturbative approach 
On 1 based on this Ansatz will improve the value of the mean energy and provide a further progress in the understanding 
of the phenomenon. 

In this note the Wannier functions are introduced by using the functions invariant under translations. This con- 
struction has some degree of arbitrariness related to the phase, which cannot be fixed a priori. Some possible natural 
^ choices are considered and numerical methods are used in order to decide which is the better one. 
^ Finally the Wannier functions are used in order to reformulate the Ansatz for the trial state. This analysis shows 

Q that there is at least another natural candidate for the trial state. By considering the particular case of filling factor 
O 1/3 numerical evidence is given that the original proposal gives the better choice. 



II. THE FORMALISM 

In this section some essential points developed in I and III are briefly reviewed. The Hamiltonian for the single 
electron is 
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(1) 



where the unit of length is 



and the unit of energy 



The complex notation is used for the vectors 



eB 

Tiijjc = — . (3) 
mc 



1 



W = Wi + iW2. 

Step down and step operators are defined 



They obey the algebra 



and moreover 



Let ifiuL be the solution of 



One gets easily 



[a, a^] = 1 [b, 6^] = 1 
[a, 6^] = [a,b]=Q 



Vn, (r) = (27™^!)-^ (^TT j '^''P V T 
It is convenient to introduce the coherent-state operator 

S{w) = exp (^^—j=(w*b — wb'^)^ — exp x r • exp {widx + W2dy) 

i.e. it is the product of a translation and a phase. It has the properties 

[S{w),H] = 

and 

S{c)S{d) = S{c + d) exp [-^-c x d • z 

(Magnetic Translation Group (MTG)). Thus they commute 

[S{c),S{d)] — Q if c X d • z = 2™ (w integer). 

Then these operators are adequate in order to impose quasi-periodic boundary conditions on a parallelogram A 
side vectors Li, L2 

S{L2)ij = e'^^V 

provided they commute 

[S'(Li), S'(L2)] = i.e. Li x L2 • z = 27r5L [gh integer). 
Then will be the degeneracy of the Landau levels. The elements of the MTG are defined by the constraint 

[S{w),S{L,)]=Q J = 1,2 

i.e. 



w= — (—71,1X2 + ^2-^1) with ni,n2 integers. 

5L 



(17) 



Let c, d be any two vectors that satisfy eq. (^_3|) and provide a commensurate tihng of the paraUelogram, i.e. 

Li = pc 

L2 = p'c + qd 

where the coefficients are integers. One has 

gL = pqu. 

A finer lattice can be introduced generated by the vectors f and g such that 

c = Kf + t'g 
d = kT + ig 

with the condition that the coefficients are integer numbers and moreover 

f X g • z = 27r 

i.e. 



(18) 
(19) 

(20) 
(21) 
(22) 



Sets of functions can be introduced which are eigenvectors of the Hamiltonian and are invariant under translations 
given by the sites of the lattice generated by c, d or by f , g. Thus solutions of 



and 



SifWn, 

can be introduced. The solutions of the above equations are given by 



(23) 



(24) 



m,n— — 00 
+ 00 

= {PIY^ exp[-i(/xm + vn)] 



m,n— — 00 



exp 



-z • (mc + nd) x r ipn^ (r + mc + nd) 



(25) 



and similarly 



+00 



C = (5l)"* E [^(/)e-'^"]"M5(g)e-^]'>„,(r) 

ra.n— — oo 

+ 00 

= {9l)-"^ E (-)""exp[-z(ma + n/3)] 



m.n— — 00 



exp 



-z • {mi + ng) X r ipnj^ (r + mf + ng). 



(26) 



The value of the parameters are fixed by the boundary conditions in eq. (|T^). See eqs. (49), (83) and (84) of I. In 
particular it is easy to find the following relations as a consistency condition for the eqs. (|23|) and (p3) 
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27rni + /i = TTKt' + na + l' f3 

27rn2 + = TT/i't + n'a + t/?. (27) 
The norm of the wave function is (see eqs. (69) and (91) of I) 

\\<f^ E / rf'r-{[5(/)]™[5(g)]'>„,(r)}>„,(r)e*(-+"«. (28) 



7n.n— — oo 



One gets 



"ff= E (-)'"V('""+"«exp(--|m/ + 71.912). (29) 



m.n— — oo 



A similar result is valid for the other set of functions 



m\'= E (-)"""e*(™^+"'')exp(-i|mc + H')- (30) 



rn.n— — oo 



The set of functions {(/i"^} form a complete set. We use capital letters to denote the normalized functions 



For computation of the matrix elements of the Coulomb part (in II and III) it was convenient to introduce a new 
set of functions {$"(^}. One considers a reference wave function ^^l^" and defines (see Appendix in III) 

^ S{w^p)<i>Zf° (32) 

where Waf3 is a standard set of vectors 

wo^p - (27r)-i [{13 - /5o) / - (a - ao) g] . (33) 

By using the composition rules in eq. ( p^ 

S{waf3)S{f) = S{f)S{wa/3) exp (-iw^^ X f • z) = S{f)S{wap) exp(-i(a - ao)) 
S{wa0)S{g) = S{g)S{wa0)exp{-iWaf3 x g • z) = S{g)S{waf3) exp{-i{(3 - (3o)) 



(34) 



it is easy to show that l>^^ satisfies eqs. (p3) and therefore differs from by a phase. 



III. WANNIER FUNCTIONS 



In this section Wannier functions ||8|-|lq| are introduced for a particle moving in a constant and homogeneous 
magnetic field. Various sets of Wannier functions can be introduced. By starting from {'I'^^} one can define (the 
index will be suppressed for convenience of notation) 

Wr,, = gl^ E e-*(™+"'^)$"'5. (35) 

a/3 

If one uses the set {'I'^f } a different Wannier functions is obtained 

Wr,s = gl^ E e-*(™+''^)$"^. (36) 

a/3 
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Yet another Wannier function is obtained if is replaced by 

S{w^,p)S{w^p,)^Zf° = exp(^(« - ao)(/3 - (37) 
All these definitions provide a set of orthonormal functions located around the lattice site 

X = rf + sg. (38) 

Moreover they have the property (see eq. (p^) 

5(r7 + s'g)W.,.(x) ^exp(^(r'f + s'g) x x • z)>V,.,s(x + r'f + s'g) 

^(-Y'''W<^r-r'Us-s')i^)- (39) 

One expects that any different choice of phases for the translation invariant functions provides Wannier functions 
with different shapes |9|. By numerical inspection one can chose the functions with better localization properties 
among the three proposal given in eqs. (^), (^) and (^7|). From Figs. |l|, ^ and || one sees clearly that the choice in 
eq. ( p5| ) (W) gives the best localization. 

The formalism developed in section || allows the introduction of other kind of Wannier functions, i.e. those 
associated to the lattice generated by the c, d vectors 17 1. 

W^")^ = ( — j ' e-*('"^+"'')*^'' (40) 

is expected to be localized around the lattice site 

X = mc + nd. (41) 

Fig. ^ shows the shape of the Wannier function Wm}n- It should be noticed that the functions {Wm.L} do not form 
a complete set of functions. 

It is interesting to elaborate further the difference between the sets {W} and {W^"^}. It can be shown (see eq. (90) 
of I) that 

fixed fj,iy 
a/3 

On the other side, {W^ ,j} evaluated on a c, d lattice site 

r = rriK + nn' 

s = mi' + UL, (43) 



gives (see eq. (^7])) 



ra -\- s(3 = m{Ka + t'/3) + nin'a + t/3) 

= m[/i — TTKt' + 27rni] + n[i' — ttk'l + 2Tm2\ (44) 



and therefore 



fixed 

Wr^s = exp(-i7r(mKi' + jik'l)) i — e-'("^+"") \ ^ \ (45) 



on the sites of the c, d lattice. 

The conclusion of this section is that one has two good candidates ({Wr,s} and {Wm.li}) which can be used for the 
construction of a trial state for the many-body problem. 
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IV. TRIAL STATE FOR FQHE 



The original Ansatz for the FQHE trial state presented in the papers II and III was formulated in terms of the wave 
function 5'^" of eqs. ( |25| ) and (|3l|). For the filling factor 1/u it is just the Slater determinant of the Ql/u functions 
with riL = 0. In the second quantization formalism where ^ is fermion field (spin is neglected), the trial state is 

Ifi) =]J /" d2a;^,t(x)*^''(x)|0). (46) 
is a unitary transformation, thus, up to a phase, the state can be written in terms of Wannier 

[ d^xV/(x)w(r,l(x)|0) (47) 



The equation 
functions 



where m, n are the sites of the lattice generated by c, d. 

From the discussion at the end of section III it follows that one can similarly construct a trial state in terms of the 
Wannier functions Wr,s given in eq. (pq) 



1^) = n / rf'a;^^(x)Wr,4x)|0) 



(48) 



Both Ansatz have as rationale the strategy of placing the N = g^/u electrons on a regular lattice. Moreover a peaked 
distribution of the electron density around the sites diminishes the energy coming from the Coulomb repulsion. 

Before examining the numerical consequences of the choices made in eqs. (^) and (p8[), it is interesting to look 
closer at the states occupied by the single electrons. In the state |0) the electron with fiiy quantum numbers is in a 
state given by (see eq. ([4^)) 



fixed fiiJ 



fixed 



(49) 



C(/3 a/3 

while in the state (see eq. (^5|)) is described by the wave function 

fixed 
a/3 



(50) 



Thus in the first case the probability amplitude for the state ^"'^ is proportional to ||'/'"'^|| and in the second is equal 
for all af3 (fixed iiiy). 

The two Ansatz given in eqs. (E^ and (Esl) can be tested in the self-consistent equation based on the Hartree-Fock 
approximation (see eq. (158) of III). Fig. ^shows that the state \ given in eq. ( |47| ) yields a better trial state (to a 
good approximation, it is unchanged by the self-consistent procedure). 
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FIG. 1. 3D plot of the modulus of a Wannier function W given in eq. (35). The lattice is regular and triangular 
(I/I = Isl ~ 2.69 in units of magnetic length). This picture is a zoom in on a domain whose size appears in Figs. 2 and 
3.. 



FIG. 2. 3D plot of the modulus of a Wannier function VV given in eq. (36). See Fig. 1 for details. 



FIG. 3. 3D plot of the modulus of a Wannier function obtained by using the representation of the MTG given by the functions 
in eq. (37). See Fig. 1 for details. 

FIG. 4. 3D plot of the modulus of a Wannier function w'"-' given in eq. (40) ((c, d) lattice, regular and triangular, u = 3 
and therefore |c| = |d| ~ 4.66). This picture is a zoom in on the domain of definition, whose size appears in Figs. 2. and 3. 

FIG. 5. Convergence patterns in the iteration procedure based on the self-consistent equation in the Hartree-Fock approxi- 
mation. "Times" refer to the Ansatz based on the functions W'"' (eq. (47)) and "circles" to the Ansatz based on the functions 
W (eq. (48)). Filling factor is 1/3. 
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